Pulling and Pushing a Cargo With a Catalytically Active Carrier 
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Catalytically active particles suspended in a liquid can move due to self-phoresis by generating 
solute gradients via chemical reactions of the solvent occurring at parts of their surface. Such 
particles can be used as carriers at the micro-scale. As a simple model for a carrier-cargo system we 
consider a catalytically active particle connected by a thin rigid rod to a catalytically inert cargo 
particle. We show that the velocity of the composite strongly depends on the relative orientation of 
the carrier-cargo link. Accordingly, there is an optimal configuration for the linkage. The subtlety 
of such carriers is underscored by the observation that a spherical particle completely covered by 
catalyst, which is motionless when isolated, acts as a carrier once attached to a cargo. 

PACS numbers: 89.20.-a, 82.70.Dd, 07.10.Cm 



Introduction. The increasing interest in the devel- 
opment of "lab on a chip" devices has led to a stringent 
need of scaling standard machinery down to micro- and 
nano-scales. This reduction in size has raised a number 
of challenging issues, such as to endow small objects 
with the capacity to perform autonomous, directional 
motion (J-Q. Towards this goal two main routes arc 
pursued. The first one consists of designing artificial 
mechanical "swimmers" by mimicking the sophisticated 
locomotion strategies of natural micro-organisms such 
as E. Coli or Spiroplasma (see, e.g., Refs. and ref- 

erences therein) . The second approach consists of trans- 
forming chemical free energy into mechanical work by 
employing phoretic mechanisms, i.e., motion induced by 
interfacial interactions @, Q ■ In the following we focus 
on this latter approach. 

Several proposals for such catalytic self-propellers 
have already been tested experimentally (see, e.g., Refs. 
P, Ulul and Refs. H Q for reviews). The un- 
derlying idea [l[ is that an asymmetric decoration of 
the surface of a particle with a catalyst, which pro- 
motes an activated reaction in the surrounding liquid 
medium, leads to a non-uniform distribution of prod- 
uct molecules around the surface of the particle. This 
non-uniform distribution gives rise toparticle motion 
through a variety of mechanisms [l|, US [HI. When the 
size of the particle is decreased towards the micron scale 
or below, viscous and surface forces start to dominate 
and inertia-based mechanisms such as bubble ejection 
propulsion become ineffective Q. (However, alternative 
propulsion mechanisms based on bubble formation can 
remain active, as for the catalytically active tubes pro- 
posed in Ref. fiof.) If the product molecules remain 
dissolved in the surrounding liquid medium, concentra- 
tion gradients develop along the surface of the parti- 
cle. It has therefore been argued @, [ll| that in such 



cases the motion of the catalyst-covered active particle 
is phoretic, i.e., it results from the interfacial interac- 
tions between the particle and the non-uniformly dis- 
tributed product molecules generated by the chemical 
reaction. Recent experimental [2, 6, la-lR 
ulation [1J, Ml , and theoretical [ill UJ1 UA 




ies of such systems have contributed to a significantly 
improved understanding of this self-induced phoresis of 
isolated active particles. 

One of the envisioned applications of such self- 
propellers is to use them as active carriers for colloidal 
transport and assembly [l(J H3] • While in classic phore- 
sis [2l], [22[ the concentration gradients are externally 
imposed and maintained, for self-phoresis the gradients 
are dynamically generated by the catalytic reaction oc- 
curring on parts of the particle surface. The interplay 
between these concentration gradients and the ensuing 
hydrodynamic flows not only powers the phoretic mo- 
tion of such active particles, but also influences the ef- 
fective interaction among each other, with nearby inert 
particles, or with bounding walls [H, [l8|- Accordingly, 
the performance of such active particles as carriers is 
expected to depend strongly on these effective interac- 
tions. 

Here we study a simple model for the transport of an 
inert particle, posing as a cargo, by an active diffusio- 
phoretic carrier 1 . The model consists of two spherical 
particles, with the surface of one of them partially cov- 
ered by a catalyst and being connected by a thin rigid 
rod. Similar to previous studies [ll|, [Jf -flpl. w e adopt 
the standard theory of diffusio-phoresis 2l[|22| and de- 
rive, within the constraints of this approach [l8j], the 
velocity of the carrier-cargo composite for two distinct 
orientations of the link [see Figs. [Tfa)-(b)]. The results 
show that the performance of the carrier, i.e., the result- 
ing velocity of the composite system, strongly depends 
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1 "Diffusio-phoresis" is restricted here to phoresis in gradients 
of electrically neutral solutes (often also denoted as "chemo- 
phoresis"). 



2 




FIG. 1: Carrier-cargo system consisting of a hard spherical carrier (top circle) of radius R with its surface Si covered 
partially (a) , (b) or completely (c) by a catalyst (hatched area) over a cap- like region. The carrier is linked via a thin rigid 
rod of length I = 2(d — R) to a hard spherical, and inert cargo with surface £2 and radius R. The product molecules (small 
hatched circles) with number density p(r, t) and of diameter a are emerging from the solvent (not shown) via a chemical 
reaction. 



on the orientation of the link, which is a paradigmatic 
consequence of the self-phoresis mechanism of the mo- 
tion. As an extreme case we study how a spherical 
particle, which is completely covered by catalyst and 
thus is unable to move on its own, can act as a carrier 
once it is attached to a cargo [see Fig. [He)]. We note 
that the emergence of motion in this latter case can be 
intuitively understood by realizing that in the limit of a 
very short linkage this configuration is equivalent with 
the so-called dimer model of an active particle [Fo| . 

Model. The system we consider is shown and char- 
acterized in Fig. [T] We assume a uniform distribution, 
with the areal density er, of the catalyst over a hemi- 
sphere of the carrier particle. The infinitesimal thin rod 
enforces that there is no relative motion between the 
carrier and the cargo. Furthermore, the rod is aligned 
such that the carrier-cargo composite possesses axial 
symmetry. The co-moving coordinate system has its z 
direction along the axis of symmetry and the xy plane 
is the midplane between the two spheres. If the carrier 
is partially covered by catalyst, there are two possible 
configurations with axial symmetry: the active area of 
the carrier points towards [Fig. QJa)] or away from [Fig. 
Gib)] the cargo. For this system the only possible mo- 
tion is a translation along the z-direction. The purpose 
of these simplifications is to sift out the main phenom- 
ena of interest, i.e., the dependence of the self-phorctic 
velocity on the internal orientation of the carrier-cargo 
system and, for the configuration shown in Fig. HTc) 
(which can be seen as a dimer model like in Ref. jl5|). 
the emergence of motion solely from linking together an 
otherwise motionless active particle with an inert cargo. 

The catalyst promotes the chemical conversion of the 
solvent surrounding the composite system into prod- 
uct molecules. Here we focus on the particular case in 



which the chemical conversion A ™ A' + B of a solvent 
molecule (A) leads to two product molecules (A 1 and 
B) only, one with properties and size being similar to 
the solvent itself (A' ps A), the other one (B) being sig- 
nificantly different 2 . Accordingly, in the following only 
this latter one is denoted as a product molecule and 
plays the role of a solute in the solvent. The product 
molecules of diameter a diffuse into the solvent charac- 
terized by a diffusion coefficient D. Thus the net result 
of the chemical conversion can be viewed approximately 
as the generation of a solute only by an ensemble of in- 
dependent sources uniformly distributed over a part of 
the carrier surface. The reaction rate vb at such a cat- 
alytic site, i.e., the number of product molecules created 
at that site per unit time t, is assumed to be indepen- 
dent of time. The number density p(r, t) of product 
molecules is considered to be so low that among them- 
selves they behave like an ideal gas. However, there is 
an interaction potential between the product molecules 
and the moving particle, which, inter alia, realizes the 
impermeability condition at the particle surface. The 
interactions between the product molecules and the sol- 
vent are accounted for effectively via the Stokes - Ein- 
stein expression D = fcsT/(37r/xa) [23[, where ks is the 
Boltzmann constant, T is the temperature, and n is the 
viscosity of the solution (i.e., solvent plus solute). 

The direct and the solvent mediated particle-solute 
interactions combined with the non-uniform distribu- 



2 It has been argued that this is approximately the case for the 
Pt catalyzed decomposition in aqueous solution of hydrogen 
peroxide (A = H2O2) into water { A' = H2O) and oxygen (B 

= o 2 ) SI El EI. 
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tion of the solute, created at the catalyst covered re- 
gion, lead to an unbalanced osmotic pressure along the 
surface of the particle which induces flow of the sol- 
vent HI [H El IH |H H and motion of the particle. 
For micromcter-size particles typical self-phoretic ve- 
locities arc of the order of /jm/s. Accordingly, both the 
Reynolds number Re ~ PsoivVR/ u (where p so lv is the 
mass density of the solvent, R the linear extension of 
the particle, and V the velocity of the translational or 
rotational motion) and the Peclet number Pe ~ VR/D 
are small [25| , Thus the hydrodynamic description ap- 
proximately reduces to the Stokes equations and the 
convection of the solute can be disregarded compared 
with its diffusive transport. 



Phoretic slip. The effective interactions between 
the particles and the self-generated, asymmetric, non- 
uniform solute number density p{r,t) around the par- 
ticles [HI, Hi], [22[ have a typical range <5, which is com- 
parable with the solute diameter a, and induce flow of 
the solution relative to the particle. In steady state 
this hydrodynamic flow in the thin surface layer ~ 6 is 
accounted for by a (phoretic) slip- velocity [2l|, [22[ , 

v s (r p ) = -6(r p )V s p(r p ) , for r p e S (= Si U E 2 ) , (1) 

where r p is a point P on the surface E = Ex U E2 of the 
particles and V s denotes the projection of the gradient 
operator onto the corresponding local tangential plane 
of the surface of the particles 3 . This serves as a bound- 
ary condition for the hydrodynamic flow in the outer 
region, i.e., outside the surface layer [H, [H,H3- The 
mobility coefficient b(r p ) is determined by the afore- 
mentioned effective interaction potential and by the hy- 
drodynamic boundary condition (stick or slip) on the 
surface of the particle, which is the inner edge of the 
surface layer [22|, [24|. For simplicity, here we take b(r p ) 
to be constant over the whole surface E. We note that 
b ^ correspond to repulsive and attractive effective 
interactions, respectively flU r22j. 
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Distribution of product molecules. In the limit 
of small Peclet numbers and neglecting any so-called 
polarization effects of the surface layer |22| > tne steady 
state distribution p(r) of product molecules in the outer 
region around the two moving particles is governed, in 
the co-moving frame, by the diffusion equation 

DVV(r) = 0, r e outer region. (2) 

This is subject to the boundary conditions (BCs) 



p(|r|->oo) = 0, 



(3a) 



3 Rigorously, Eq. {T} should be interpreted as a condition on 
the outer edge of the surface layer. For the outer problem 
this can be replaced by E because p varies over length scales 
which are much larger than 8. 



1, re catalyst 
, otherwise 

(3b) 

Equation (|3b[) , where vb <t corresponds to the reaction 
rate per unit area and n denotes the outward direc- 
tion normal to the surface, describes how the catalytic 
reaction at r £ E translates into a source of product 
particles B [UGl]. 

Hydrodynamic flow outside the surface layer. 
Because there are no forces acting on the solution be- 
yond the surface layer, based on the assumption of low 
Reynolds numbers the hydrodynamic flow field u in the 
outer region is obtained as the solution of force free and 
incompressible Stokes equations: 



v • n = 0, v • u = 0. 



(4) 



II := —pi + pS is the corresponding pressure tensor, 
where p is the hydrostatic pressure and S is the shear 
stress tensor, i.e., S a p = du a /dxp + dup/dx a . In the 
comoving reference frame, attached to the carrier-cargo 
composite which is moving with V = Ve z relative to the 
solution, which is quiescent far away, these equations 
are subject to the BCs of a flow velocity —V far away 
(|r| — s- 00) from the composite and v s (r p ) at the outer 
edge of the surface layer. This corresponds to sticking 
at the surface of the particles plus a slip velocity v s at 
the outer edge of the surface layer: 



|r|- 



V. 



(5) 



Equation ([5]) reveals that the cargo plays an active role 
(v s 7^ on £2 C E) despite being catalytically inert. 

Phoretic velocity. After computing the hydrody- 
namic flow in the outer region, which depends para- 
metrically on the translation of the two particles via 
the BCs in Eq. ([5]), the phoretic velocity V is deter- 
mined by the condition that the motion of the system 
composed of the particles plus their surface layers is 
force free [22| . 

Equations @ and (J3J) as well as Eqs. ^ and are 
most conveniently solved in terms of bispherical coor- 
dinates (£ e k, < 77 < 7T, < 4>< 2tt) Hi-ii] 

{x , y , z} = >r{sin?7 cos <fi , sin 77 sin 4> , sinh£}/C , (6) 



where C = cosh £— cos 77; the corresponding scale factors 

1 H = \8r/dt\ = [{dx/di) 2 + {dy/d0 2 + (dz/di) 2 ] 1/2 , 
hr) = \dr/dn\, and h# = \dr/dc/>\ are given by 

{h,£ , h n , h^} = x{l , 1 , sin77}/C. 



(7) 



The choice x = v 'd 2 — R 2 ensures that the family of 
spheres £ = const includes the ones, £ = ±£0 := 
± arccosh(<i/i?), corresponding to the surfaces Ei j2 of 
the carrier and the cargo, respectively. Noting that the 
intersection of an rj iso-surface with a £ iso-surface is 
a circle parallel to the xy plane and that the equato- 
rial circle parallel to the xy plane on Ei corresponds 
to 770 = arcctg(i?/ \fd 2 — R 2 ), the area covered by the 
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catalyst is parametrized by £ = £o , < <f> < 2w, 
and tjq < r\ < ir for the configuration in Fig. HJa), 
< r\ < ?7o for the one in Fig. [TJb), and < r\ < ir for 
the one in Fig. [He). 

In terms of bispherical coordinates the solution p(£, 77) 
of Eq. © , which does not depend on </> due to the axial 
symmetry of the system, which is finite at rj = 0, 7r (i.e., 
on the z axis), and which satisfies the BC in Eq. (|5a|) . 
can be written as [2a, u3 



p(£,oj := cos??) = p (cosh£ - cj) 1/2 ^ P n (u) 

n>0 

x [A„ sinh ((n + l/2)£) + P« cosh ((n + 1/2)0] . (8) 

where P n denotes the Lcgendrc polynomial of or- 
der n and po := Rvbo/D is a density scale 
chosen such that the system specific parameters 
are factored out from the boundary conditions 
[Eqs. Noting that n = (-,+)e 4 = 

(+, — ) {sinh £ sin r\ cos <p , sinh £ sin 77 sin </> , cosh £ cos 77 — 
1}/C on £^2 and expanding the right hand side of Eq. 
(|3b[) in a series of Legendre polynomials, the BCs in 
Eq. (|3b[) determine the coefficients {A n ,B n } as the 
solution of a system of linear equations. This deter 
mines p(£, 77) and the slip velocity follows from v 
-[(b/h TI )d v p(£,ri)]£=±£ with e r) = {(cosh£cos?7 - 
1) cos 4> , (cosh £ cos 77 — 1) sm(j> , — sinh£ sin^j/C. 

Due to the linearity of the Stokes equations the hy- 
drodynamic flow field u = uj + ujj [Eqs. @ and 
([5])] is the superposition of the one corresponding to 
Uf(|r| — > 00) = —V with stick BCs on S and the 
one corresponding to a quiescent flow far away, i.e., 
ur/(|r| —> 00 = with slip BCs v s on S. The so- 
lution for the first problem is known [28[ and cor- 
responds to a hydrodynamic force on the composite 
F/ = -127r/ii?AV, where 



A 



n(n + 1) 



sinh £n , 

^ ^ (2n- l)(2n + 3) 

n>l 7 

4sinh 2 (n+ l/2)£ - (2n + l) 2 sinh 2 £ 
2 sinh(2n + 1)£ + (2n + 1) sinh 2£ 



(.9) 



The second solution u// for the flow field is most conve- 
niently expressed in terms of the Stokes stream function 
\Jj for axisymmetric flows [which in cylindrical coordi- 
nates {f = y^x 2 + y 2 , 0, z} is defined by (uf,u z ) := 
(—d z ^,df-^)/r\. In terms of bispherical coordinates 
*$>(£, uj := cos 77) is given by (28| 



n>0 



(cosh£ — uj) 3 / 2 



(10) 



mally denoted as the Gegenbauer polynomial of degree 
n and order —1/2 [30j . and the characteristic velocity 
Vo = bpo/R is chosen such that the system dependent 
parameters are factored out from the boundary condi- 
tions in Eq. flSJ. Because there are no flow sources 
or sinks along the z-axis, we require the stream func- 
tion to vanish there (25|, i.e., w = ±1) = 0. Since 
Cm 1 (uj — ±1) = for m > 2, but as a linear func- 
tion c\ 1 ^ 2 \uj) cannot be zero for both uj = ±1, this 
constraint implies that the term n = is removed from 
the series representation in Eq. (JTUJ). 

In terms of the stream function, the boundary 
condition of a vanishing flow velocity along the di- 
rection normal to the surfaces Si. 2 (correspond- 
ing to impenetrability) leads to the requirement 
[9\I>(£, u = cosr])/dr] = 0] ± ^ Q = (25[, which can be 
shown to be equivalent to the simpler form ^(±£0, oj) = 
[HI, [3(|. The flow velocity along the slip di- 
rection 7] is given in terms of the stream function 
by u v = —C 2 (x 2 sin?;) -1 ^^ [2f| which, according 
to Eq. ([5]), on the surfaces Si i2 should be equal 
to the phoretic slip v s . Since for n =/= m and 
n,m > 2 the polynomials Cn,m 2 ^ 



- 1/2) Mc& 



1/2) 



are orthogonal, 
(uj) = 0, the co- 



efficients {a n ,...,d n } can be determined from these 
boundary conditions on ^ in terms of {A n ,B n } by 
solving a system of linear equations 4 . The hydrody- 
namic force acting in the positive z-direction on the 
composite due to this second flow u/j is given by 
F n = (2 5 / 2 irpV R 2 /K)Y /n > 1 (a n +c n ) @. 

Requiring that the total force F/ + F/j acting on the 
composite vanishes [2l], [22j thus leads to the following 
expression for the velocity: 



V/V 



V2 



3X v / (d/R) 2 - 



= *^2( a n + C n ) ■ 



(12) 



n>l 



Discussion. In Fig. [2] we show the velocity of 
the carrier-cargo system as a function of the scaled 
length l/R = 2 (d/R — 1) of the linkage in the range 
0.02 < l/R < 18 for the three configurations considered 
in Fig. [T] in the case of repulsive effective interactions 
between the product molecules and the material form- 
ing the carrier-cargo system, i.e., for b < and Vb < 0. 
(The calculations can easily be extended to l/R < 0.02, 
but: (i) the implicit assumptions of a homogeneous sol- 
vent and point-like product molecules are expected to 
break down for very small separations, which renders 
Eq. ([T2l to be inapplicable for too small values of I /R, 
and (ii) the limiting configuration l/R — correspond- 



where 



W„(0 = a n cosh ((n - l/2)£) + b n sinh ((n - 1/2)0 
+c„ cosh ((n + 3/2)0 + d n sinh ((n + 3/2)0 , (11) 

C ( - 1/2) (uj) = [P n _ 2 (w) - Pn(w)]/(2n - 1) H is for- 



4 In practice, the infinite system has to be truncated at a certain 
order n = N max . For our system it turns out that both An 
and B n vanish rapidly upon increasing n. In all cases stud- 
ied, at N max = 100 these coefficients are smaller than 10 -20 , 
compared with values of the order of unity for n = 0. This be- 
havior is sufficient for the convergence of the velocity expression 
in Eq. JT3J. 
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FIG. 2: The velocity V along the ^-direction (indicated in the insets) in units of \Vo\ as a function of the scaled length 
0.02 < l/R < 18 of the linkage for the three configurations shown in the corresponding insets and for the case that the 
effective interaction between the product molecules and the two particles is repulsive (b < 0). The dashed lines in (a) and 
(b) correspond to the values ±|Vo|/8, respectively, which the velocity of the cargo-carrier composite seems to attain as the 
cargo is linked at a larger distance from the carrier. 



ing to carrier and cargo being in contact merely reduces 
to the case of self-diffusiophorcsis of a single particle 
with a complex shape.) 

For the choice b < the carrier alone moves against 
the gradient of product molecule along its surface. 
Therefore the configurations in Figs. Q] (a) and (b) 
correspond to a pulling and a pushing carrier, respec- 
tively. Focusing first on the configurations in Figs. [2] 
(a) and (b), we note that the pusher or puller charac- 
ter is maintained at all separations I, i.e., the velocity 
of the composite is always along the positive and neg- 
ative z direction in the cases (a) and (b), respectively 
In all cases, for l/R < 4 the velocity shows a signifi- 
cant dependence on the separation I, which clearly in- 
dicates that the cargo plays actually an active role in 
the resulting motion of the composite. In both cases 
at large separations / the velocity approaches limiting 
values seemingly equal to ±|Vb|/8, respectively, which 
is 1/2 of the velocity at which the active carrier half- 
covered by catalyst would move if isolated [l^. How- 
ever, the pusher slows down upon increasing I whereas 
the puller speeds up. Therefore, for any length I of the 
linkage the pusher configuration exhibits a faster mo- 
tion than the puller one. These findings can be ratio- 
nalized intuitively as follows. For example, in the case 
(a) the generation of product molecules by the carrier 
induces also a concentration gradient around the cargo 
surface pointing in the direction opposite to the one on 
the carrier surface. Consequently, the cargo behaves 
also like an active particle and induces a hydrodynamic 
flow around its surface in the direction opposite to the 
one around the carrier, which thus leads to a smaller net 
hydrodynamic flow and accordingly to a slower motion 
of the composite. 

Any other orientation of the carrier-cargo link breaks 
the axial symmetry, inducing a rotation of the compos- 
ite in addition to its translation; this leads to a reduc- 
tion in the translational velocity. Therefore the pusher 
configuration [Fig. [T] (b)] is the optimal (i.e., maxi- 
mal velocity) configuration if the effective interaction 
between the product molecules and the carrier-cargo 
material is repulsive. On the other hand, if this ef- 
fective interaction is attractive (b > 0) the same intu- 



itive arguments as above lead to the conclusion that the 
puller and pusher characters are interchanged between 
the configurations shown in Fig. [lja) and (b), while the 
/-dependence of the velocity changes signs. Therefore, 
for attractive effective interactions the puller configura- 
tion [Fig. [T] (b)] is the optimal one. 

The motion of active particles is known to be affected 
by thermal fluctuations, and thus its directionality only 
holds over time scales smaller than the characteristic 
rotational diffusion time, cross ing over to Brownian mo- 
tion at longer time scales 0, [3l[ . Inspecting the cases 
of the carrier-cargo configurations shown in Figs. [1] (a) 
and (b), one realizes that in (a) the cargo has a stronger 
influence than in (b) on the stability of the compos- 
ite with respect to such fluctuations of the direction of 
motion. For b < 0, in the case of the puller [Fig. [1] 
(a)] a sudden change in direction producing a tilt of 
the symmetry axis towards the left relative to the Oz 
axis breaks the axial symmetry of the product concen- 
tration profile, and larger concentrations occur to the 
right of the symmetry axis. Consequently, right after 
the change in direction the cargo is driven through a 
region where the concentration of the product particles 
is larger on its right half (relative to the cargo-carrier 
axis) and thus induces a torque on the cargo and a ro- 
tation of the symmetry axis back towards the Oz direc- 
tion. Therefore, in the puller configuration the cargo 
tends to stabilize the directionality of motion against 
thermal fluctuations. On the other hand, in the pusher 
configuration [Fig. [U(b)] the cargo is further away from 
the source of product particles. Therefore it is basically 
insensitive to such changes in the symmetry of the prod- 
uct concentration distribution as it moves away from it 
and consequently no significant restoring torques occur. 

Configuration (c) reveals the rather peculiar conse- 
quence of self-phoretic propulsion that a carrier, which 
in isolation would be motionless, is activated by be- 
ing linked to an otherwise inert cargo. This occurs be- 
cause the linkage provides the anisotropy needed for 
self-phoresis to become operational 0, [TiJ. For short 
links, the velocity of such a composite is significant. 
According to Fig. [2 for l/R < 1, |V|/|Vb| is approx- 
imately 1/3 of the largest possible velocity |Vb|/4 of 
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an isolated spherical active carrier. For repulsive (at- 
tractive) effective interactions the inert cargo turns the 
neutral carrier into a pusher (puller). We note that for 
lengths I /R <C 1 the carrier-cargo system in configura- 
tion (c) becomes similar to the geometry of the so-called 
dimcr model of an active colloid which, for an A — > B 
catalytic reaction, has been proposed and investigated 
with Molecular Dynamics simulations [l5j . 

The diffusiophoresis mechanism and the geometry of 
the model have been chosen with the intention to pro- 
vide clear examples which reveal and underscore the 
complexity of such active carrier - cargo systems. We 
note that the main influence of the inert cargo stems 
from the fact that the nearby chemically active col- 
loid gives rise to concentration gradients, and thus to 
phoretic slip velocity, along its surface, too. This is very 
different from the cases discussed in Ref. [32| , where hy- 
drodynamic interactions between dimer swimmers are 
" rectifying" the otherwise reciprocal movement of each 
of two dimers, and thus are leading to the emergence of 
collective and relative motion. We anticipate a rather 
rich behavior to emerge as various model constraints are 



relaxed, such as: (i) Different effective interactions be- 
tween the product molecules and the cargo and carrier 
material, i.e., different mobility factors b (of potentially 
opposite sign) on the carrier and cargo surfaces, will 
lead to a complex dependence of the velocity on the 
orientation of the linkage. (For example, in the case 
shown in Fig. [21(a) with bx,i of opposite sign the veloc- 
ity can change sign as | £»2 1 is increased.) (ii) In the case 
of charged active particles and charged reaction prod- 
ucts, the motion will be determined by the interplay of 
self-diffusiophoresis and self-electrophoresis, which am- 
plifies the roles played by the orientation of the linkage 
and by the surface properties of the cargo and carrier, 
(iii) A convex or concave shape of the surface of the 
cargo facing the carrier side leads to a possible ampli- 
fication or reduction, respectively, of the concentration 
gradients around the carrier and the cargo. This list can 
be extended, but it already now shows that employing 
active particles as carriers allows an exceptional flexibil- 
ity in the design of cargo-carrier systems, which should 
be very beneficial for potential applications. 
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